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We exactly solve the ferromagnetic spin-1/2 Ising model on the Bethe lattice in the presence of 
an external magnetic field by means of the equations of motion method within the Green's function 
formalism. In particular, such an approach is applied to an isomorphic model of localized Fermi 
f*"^ ■ particles interacting via an intersite Coulomb interaction. A complete set of eigenoperators is found 

together with the corresponding eigenvalues. The Green's functions and the correlation functions 
are written in terms of a finite set of parameters to be self-consistently determined. A procedure 
is developed, that allows us to exactly fix the unknown parameters in the case of a Bethe lattice 
with any coordination number z. Non-local correlation functions up to four points are also provided 
together with a study of the relevant thermodynamic quantities. 
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<D I I. INTRODUCTION 

Recently, it has been shown [l[ that a system built up of q species of Fermi particles, localized on the sites of 
a Bravais lattice and subjected to finite-range interactions, is exactly solvable in any dimension. Exactly solvable 
means that it is always possible to find a complete set of eigenvalues and eigenoperators of the Hamiltonian, which 
-4_J ■ close the hierarchy of the equations of motion. In such a way, exact expressions for the relevant Green's functions 
and correlation functions can be derived. These expressions areiust formal because they depend on a finite set of 
parameters to be self-consistently determined. In Refs. @, d, [j] it has been shown how it is possible to fix such 
parameters exactly by means of algebra constraints in the case of one dimension and q = 1, q = 2 and q = 3, 
y , respectively. In this way, complete and exact solutions of these systems have been obtained. 

A system of q species of Fermi particles has been shown [l[ to be isomorphic to a spin-| Ising-like model in the 
[ presence of an external magnetic field, so opening a different route to the study of spin systems, which can be very 
l— ~~ '■ difficult in two and three dimensions when attached by the transfer matrix method. Furthermore, this approach can 
shed new light on how to get an exact solution for these systems in higher dimensions in the presence of an external 
magnetic field as it is always possible [l| to find an exact expression for the corresponding Green's functions and 
' correlation functions. The exact knowledge of the eigenenergies of the system can give information on the energy 
Q\ , scales ruling the physical behavior and the response of the system and can find an application as unbiased check for 
^sO • the approximate studies present in the literature. Within our approach the problem is that correlation and Green's 
functions depend on a finite set of unknown parameters to be self-consistently determined. A complete exact solution 
\Q • of the system is obtained only when such parameters are known. These parameters cannot be determined by means 
of the dynamics and are fixed by choosing the representation where the field operators are realized. In particular 
C . [cfr. Section 2.4 in Ref. 0], they can be fixed by appropriate self-consistent equations which are the manifestation 
of symmetries of the model, algebraic properties of the field operators, boundary conditions (i.e. properties of the 
^ , underlying lattice, phase of the system according to the values of the external thermodynamical parameters). It has 
been shown how to fix exactly such parameters by means of algebra constraints in the one dimensional (ID) case 
0, [H, H[ and we are now working on the possibility to introduce new algebra constraints and topological relations in 
order to fix the self-consistent parameters in higher dimensions. The first step to realize such an ambitious program 
appears to be the application of our technique to a Bethe lattice of any coordination number z. In this article we 
apply our formulation, the Composite Operator Method (COM) [H, Q , to the study of a system of q = 1 species of 
Fermi particles interacting with an intersite Coulomb interaction on the Bethe lattice with coordination number z. 
Such a model is shown to be isomorphic to the spin-i Ising model on the Bethe lattice in the presence of an external 
magnetic field. 

Bethe Q and Bethe-like lattices Q have been widely studied in solid state and statistical physics because they 
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represent the underlying lattices on which many systems can be exactly solved [§, E3, E2, EH EH, E3, El, Ef| [13, 
El HI E3- Following the line of reasoning that refers to the mean field theory of magnetism, the Bethe-Peierls 
approximation was developed in 1935/36 in order to describe crystalline alloys or Ising models. It takes exactly into 
account the interaction of a given spin with its nearest neighbors but introduces a mean field in order to express 
the interactions between such neighbors and all the other spins in the lattice. Later, it was pointed out [231 ] that 
such an approximation becomes exact on the Bethe lattice. In particular, it was shown that the partition function 
of the ferromagnetic spin-^ Ising model on the Bethe lattice with any coordination number z is equivalent to that 
in the Bethe approximation (23[. Then, the equivalence of the exact solution of the Bethe lattice spin-i Ising 
model to the Bethe-Peierls approximation was shown also in the antiferromagnetic case [HI]. Summarizing, there 
are two special properties that make Bethe lattices particularly suited for theoretical investigations: the self-similar 
structure which may lead to recursive solutions and the absence of closed loops which restricts interference effects of 
quantum-mechanical particles in the case of nearest-neighbor coupling. Furthermore, Bethe and Bethe-like lattices 
have attracted a lot of interest because they usually reflect essential features of systems even when conventional 
mean- field theories fail [13] • The reason is that such lattices are capable to take into account correlations which are 
usually lost in conventional mean-field calculations. The spin-^ Ising model on the Bethe lattice can be exactly solved 
by means of the transfer matrix technique [25| which reduces the solution to an eigenvalue problem of the second 
order and all the relevant thermodynamic quantities such as the magnetization and the susceptibility, and so on, can 
be calculated by means of recursion relations, thanks to the nested structure of the underlying lattice. The same 
technique has been recently employed in the exact calculation of the spin-spin correlation functions (S (0) S (n)} for 
any temperatures T and external field h [26l . [27j . Exact expressions for the free energy and the magnetization of 
a spin-i Ising model on a two-layer Bethe lattice in the presence of magnetic fields different in the two layers have 
been also obtained together with a study of the whole phase diagram by means of an iteration technique [28]. Also, 
it should be mentioned a large activity in the framework of the athermal random- field Ising model (RFIM), where 
analytical results have been obtained on Bethe lattices [1^, [3(| ■ 

In this paper we exactly solve the ferromagnetic spin-^ Ising model on the Bethe lattice with any coordination 
number z in the presence of an external magnetic field within the COM approach [1, Q. All the Green's functions and 
correlation functions are obtained together with the behavior of the relevant thermodynamic properties. Two-point 
(S (0) S (j)) and three-point (S (0) S (j) S (fc)) spin-spin correlation functions are also provided together with non 
local correlation functions of higher order. The manuscript is organized as follows. In Section 2, we give the general 
Hamiltonian of the spin-i Ising model on the Bethe lattice and the mapping onto a model of Fermi particles with 
intersite Coulomb interactions. In Section 3, we present the general solution in terms of eigenvalues and eigenvectors. 
In Section 4, we show how to close the system of self-consistent equations and find the unknown parameters in order 
to compute the correlation and the Green's functions. In Section 5, we compute the local correlation functions and in 
Section 6, the non-local ones. In Section 7, we study all the relevant thermodynamic quantities, such as magnetization, 
susceptibility, internal energy, specific heat and entropy as functions of the temperature and the external magnetic 
field, specializing the general formulas to the case of a Bethe lattice with coordination number z — 3 and z = 4. 
Finally, some concluding remarks and outlooks of our work are given. Some technical Appendices follow. 

As a final remark, we would like to stress that the motivation of this work is to show that the formalism of Green's 
functions and equations of motion is a convenient technique to study spin systems. Most of the techniques used in 
the literature for the study of these systems are based on the transfer matrix method. This latter formalism is a 
very powerful technique and has been largely applied with success to a huge number of models. After the brilliant 
solution by Onsager [31| for the two-dimensional spin-i Ising model in zero field, many other two-dimensional (2D) 
models, such as the dimer problem, six- vertex, eight-vertex (see Baxter's book), have been solved by making use of 
the transfer matrix method [see Baxter's book [25| for a comprehensive list of references]. However, it should be 
noticed that this method is very transparent and convenient for the case of one dimension, but becomes complicate 
for higher dimensions. In spite of the tremendous work done [among the most recent results, the derivation of the 
order parameter of the chiral Potts model by Baxter [32| has to be mentioned], many problems remain unsolved. The 
exact partition function in a finite magnetic field is still unknown. No exact results have been obtained for the three- 
dimensional model. By using the equation of motion formalism, we have constructed a general method to study Ising 
spin systems [l| . Such a method is general, in the sense that it has been formulated for any dimension of the system. 
We can exactly calculate a complete set of eigenoperators and eigenvalues of the Hamiltonian, and consequently to 
derive analytical expressions for the correlation functions. In order that this scheme of calculation could be used in 
practice, it is necessary to calculate a set of unknown parameters. The number of unknown parameters depends not 
only on the dimensions of the system, but also on the dimension of the spin; for a Ising spin-g/2 system on a lattice of 
coordination number z, the number of unknown parameters is 2qz [cfr. Ref. [lj]. Our previous studies muni show 
that it is possible to find the necessary self-consistent equations by using not only properties of the lattice, but also 
symmetry and algebraic properties of the field operators. In the last two years we have been performing a systematic 
study of this last point. We started by considering the simplest problem of spin- 1/2 on a linear chain 0. Then, we 
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considered the case of spin-1 Q and spin-3/2 Q, always for ID systems. The extension to spin higher than 1/2 is 
not immediate, but requires the introduction of higher composite fields (projection operators). After this study of 
ID systems, we decided to consider more complicate lattices, by considering the Bethe lattice. This lattice has the 
same topology of ID because the absence of closed loops, but the analysis requires a dependence on the coordination 
number. For the Bcthc lattice we have shown that the problem can be completely solved; we have shown that all the 
known results existing in literature can be reproduced. Furthermore, we have obtained new results, not previously 
obtained, as for the case of three-point correlation functions. The next step we have in program is the study of the 
spin-1/2 for the 2D lattice. This step is a very hard task; the properties of the lattice are different and the introduction 
of new concepts for writing down the self-consistent equations for the unknown parameters will be necessary. 



II. THE MODEL 



Let us consider the spin-^ Ising model with nearest-neighbor interactions, in presence of an uniform external 
magnetic field h, on a Cayley tree with coordination number z. The Hamiltonian can be written as: 

z 

H = -hS (0) + H{P) (!) 
P =i 

where S (0) is the spin operator at the central site (0). The spin variables S take only two values: S = ±1. flw is 
the Hamiltonian of the p-th sub-tree rooted at the site (0) and can be written as: 

Z-l 

fl-(f) = -hS (p) - JS (0) S (p) + H(p,m) (2) 

m— 1 

where (p), (p — 1, z) are the nearest neighbors of (0), also termed the first shell. In turn ij(P' m ) describes the m-th 
sub-tree rooted at the site (p). The process may be continued until we eventually reach the boundary sites, described 
by the Hamiltonian: 

H^-^ = -hS(p 1 ,...,p r )-JS(p 1 ,...,p r . 1 )S( Pl ,...,p r ), (3) 

where (pi,p2, ...,Pr) [pi = 1 ; z; P2,P3, ... = 1, z — 1] are the boundary points belonging to the r-th shell. In what 
follows we focus only on the sites deep in the interior of the tree, so ignoring the boundary, i. e. we concentrate on 
the Bethe lattice. Let us now consider the transformation: 

S (i) = 2n (i) - 1 (4) 

where i is a generic site of the lattice, 

n (i) = c f (i) c (i) (5) 

is the density operator for a spinless fermionic field, c (i) and (i) being the annihilation and creation operators 
satisfying the canonical anti-commutation relations: 



{c(i,f),ct(j,t)} =dy 

{c(i J *),c(j,t)} = {ct (M),ct(j,i)} = 



(6) 



In this way a mapping is established between the spin-i Ising model and a model of Fermi particles with intersite 
Coulomb interactions on the Bethe lattice, where the correspondence between the Ising and the fermionic variables 
is: 

S = 1 => n = 1 (J s 
S = -l => n = 0" 1 ' 

The Ising Hamiltonian, eqs. CO)-©, with the transformation ((4]), takes the form: 



H = E + 2(zJ-h)n (0) + £* =1 ftW 
fiW = 2 (zJ — h)n (p) - 4 Jn (0) n (p) + Y^li # (p,m) 



(8) 

jj(pi,-,Pr) =2{zJ-K)n (p 1; ...,p r ) - AJn (p 1; ...,p r _i) n (pi, ...,p r ) 
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where the constant term Eq is defined as: 

E = h + z {h - J) ( z - !) P_1 = h + z(h — J) [Z . (9) 

P =i z 

We immediately recognize the chemical potential \x — 2 (h — zJ) and the potential strength V = —4 J in a fermionic 
language. Also here we ignore the boundary sites and reduce to the Bethe lattice. Such an Hamiltonian enjoys the 
particle-hole symmetry, that is, it turns out to be invariant under the transformation n — > 1 — n, which in the spin 
language corresponds to the spin-inversion symmetry S — > —5, h — > —h; in particular the chemical potential as a 
function of n scales as 

fj,(l-n)=zV-fjL{n). (10) 

This scaling law implies that the magnetization vanishes in zero external magnetic field. However, as it will be 
shown in Section 7 the Hamiltonian (Q]) and/or ([5]) admits also solutions exhibiting a spontaneous breakdown of the 
particle-hole symmetry; that is a magnetization different from zero in absence of magnetic field. 
We see that the density operator satisfies the equation of motion: 

i^-n(i) = [n(i),H} = 0, (11) 

so that standard methods based on the use of equations of motion and Green's function formalism are not immediately 
applicable in terms of this operator. The relevant equation of motion to be considered is: 

d 

i—c{i) = -fic(i)-4:zJc{i)n a {i) (12) 

where 

n a (i) = ± (13) 
P =i 

(i,p) being the nearest neighbors of the site i. 

In the next Section, we will show in detail how to deal with such an issue and build up the formalism. We shall 
put the attention to the fermionic system and will solve the Hamiltonian |8]) by using the formalism of the equation 
of motion and Green's function method [l], [3, S EL IE @| ■ 

III. COMPOSITE OPERATORS AND EQUATIONS OF MOTION 

In this Section, we exactly solve the Hamiltonian {8} starting from the identification of a suitable operatorial basis 
[H, @. In order to pursue this task, we focus on the central site (0), even though we could have chosen any other site 
thanks to the symmetry of the Bethe lattice. Let us consider the following series of composite field operators 

^(0)^(0)^(0)]", fc = l,2,... (14) 

where, according to the definition (fl"3|) n a (0) = | Ylp=i 71 (p)i (p) being the first neighbors of the site (0). By using 
(fTTj) and (jT3J) it is easy to see that these operators satisfy the hierarchy of equations of motion: 

(0) = bPk (0) , H] = -fiiPt (0) - 4z J^ fe+1 (0) . (15) 

at 

However, we observe that the number operator n (i) satisfies the following algebra: 

[n = [ct (<) c (<)] k = n (*) , k = 1, 2, ... . (16) 



of the fields [n a (i)] k 



As shown in Appendices A and B, this algebraic property allows us to establish the following fundamental property 

z 

KWf^^Kfi)]" 1 (17) 
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where the coefficients A$ are rational numbers which can be calculated according to the scheme given in Appendix 
B. Therefore, for k = z + 1 the hierarchy of equations of motion (|15[) closes as the additional operator (0) = 
c(0) [n a (0)] z+1 can be rewritten in terms of the previous z + 1 elements of (fl4|) through the relation (fTTj) . We are 
thus able to derive a closed set of eigenoperators of the Hamiltonian by defining the following composite operator: 



V>(0) 



/ V-i(o) \ 




( c(0) \ 


i>i (o) 




c(0) n a (0) 


V 4>*+i (o) J 




[c(0) [n a (0)f J 



(18) 



which satisfies the equation of motion: 



i-V(0) = [^ (0),i/]=eV (0) 



(19) 



where the (z + l)x(z+l) energy matrix e is defined as: 











V o 



-4zJ 





-4zJ 






-AzJA 



(«H-1) 







2-2 



-4zJ 



-4zJA 





-4zJ 



(*+i) 

z-1 



(20) 



The eigenvalues of the energy matrix have the expressions 

E„ = -fx - 4(n - 1) J, n = 1, 2, ar + 1 



(21) 



At this stage we can say that we have formally, but exactly, solved Hamiltonian J5]) or its spin counterpart (U)-© 
as we have found for them a complete set of eigenoperators and eigenvalues for any coordination number z of the 
underlying Bethe lattice. The solution is formal as we have to compute still the correlation functions. 
In order to do this, let us now define the thermal retarded Green's function 



where (...) 
transform: 



G R (t - t') = (R (0, t) (0, f )] > = e (t - H) {{^ (0, t) , (0, t')}) (22) 
denotes the quantum-statistical average over the grand canonical ensamble. By introducing the Fourier 



G R (t - t') = 



(27T) 



+00 



and by means of the Heisenberg equation (| 19(1 we get the equation: 

[u - e] G R (to) = I 

where / is the normalization matrix, defined as: 

/=({^(0,i),^ (0,i)}) 

The solution of Eq. pi]) is (H. 1^ 



LO - E n 



(23) 

(24) 
(25) 

(26) 



where are the spectral density matrices, to be calculated through the formula 0,0): 

2 + 1 



<cb, 



(27) 
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where f2 is the (z + 1) x (z + 1) matrix whose columns are the eigenvectors of the matrix e. 
The matrix f2 has the expression 



fl 



n = 



z- 1 (§)' 
(§)' 



z 2 



z> (f) 2 • 

*i (f) 1 • 

^0 z° (f)° • 
in general, the matrix element fi P) fc has the expression: 



(a) 



2-2 



z-2 



z-l 

z 
2-1 

z 

2-1 



1 
1 



(28) 



2+1 — p 



k = l,p = 1 
* = 1,P^ 1 



(29) 



By means of the definition 
has the expression 



and of the recursion rule Q17p , the normalization matrix can be easily calculated and 



I = 



( 1 


ha 


h,3 ■ ■ 


■ h.z-l 


h,, 




h.2 


h,a 


h,4 ■ ■ 


■ h.z 


h,z+i 


h,z+l 


h,3 


h,i 


h,s ■ ■ 


h,z+i 


^2,2 + 1 


h,z+i 


h,z-i 


h,z 


h,z+i ■ ■ 


' ^2-3,2 + 1 


^2-2,2 + 1 


Iz-l,z + 


h,z 


h, 2+1 


h,z+i ■ ■ 


' -^2-2,2+1 


^2-1,2 + 1 


Iz,z+1 


\ h, 2+1 


-^2,2 + 1 


h,z+i ■ ■ 


' ^2-1,2+1 


Iz,x+1 


Iz+l,z+ 


2,...,z + 


1) are 


expressed 


as 







(30) 



I p , z+1 = A%+'-»I lim+ i. 

Ill 1 

Therefore we need to know only the elements = 2, z + 1) which are given by 



K(o)] 



fc-i 



Jk-l) 



(31) 



(32) 



Then, the spectral density matrices a^ n ' can be easily calculated by means of Eq. (|27|) once we keep in mind that, 
according to the structure of the normalization matrix J, Eqs. (|30|) and (|31|) . there exist only z + 1 independent 
matrix elements crj™2 f° r each of the z + 1 matrices tr*-™- 1 while all the others can be obtained as linear combinations 
of these latter according to (|3"Tj) . As a result we get: 



(33) 



where E n are functions of the elements Ji^ with fc = 1, ...,z + 1 and are numerical matrices. In particular we 
have the following expressions: 



s i - Efcti fi i,/Ui,fc 



and 



rf> = (i,o,o,. ..,o,o,o) k = i,...,z + i 



l,fe 



{~T~) k \ n = 2,...,*+l 



(34) 



(35) 
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The correlation function 



C(t-t')=(i>(0,t)^ (0,i')> = 7^7 

can be computed starting from Eq. ff^5|) and recalling the relation: 



+00 



rfwe"" ('"*') C(u) 



C(u) 



Then we get: 



1 + tanh 



2+1 



(3uj 
~2~ 



C{uj) = nJ2° in)T nS[u-En\ 

n=l 

C (t-0 = i^e-'^^-OtrWTn 



where 



1 + tanh 



V 2 



(36) 
(37) 

(38) 
(39) 

(40) 



Eqs. (f!?6")) and (|38p are the exact solution of the problem. Such a solution is only formal because the complete 
knowledge of the retarded and correlation functions is not fully achieved owing to the presence of the unknown 
static correlation functions — ([n a (0)] m ) (m = 1, ...,z) appearing in the normalization matrix I due to the 

non-canonical algebra satisfied by the composite field operator ip (0). Such unknown parameters will be calculated 
according to the self-consistent scheme given in the following Section. 



IV. SELF-CONSISTENCY 



As we have shown in the previous Section, a complete solution of the model requires the knowledge of the correlators 
(0)]™) . In order to compute these quantities, let us write the Hamiltonian (jSJ) as the sum of two commuting 



K (m) = ^ n a 

terms: 

H = Hn + Hj 



Hi = -4Jn(0)J2 n (p) ■ ( 41 ) 



Because [Hq, Hi] = 0, for any operator O we can write its average as 



Tr{e-P R } Tr{e-P H 'e~P H °} (e 
where (...) is the trace with respect to the reduced Hamiltonian Hq 

_ Tr{...e-P H °} 



/o 



(43) 



Let us now consider the correlation functions Ci.k = (0) (0) [n a (0)] k y , k = 1, z + 1. By means of Eq. (f42|) 
we can derive the following relation: 

= V 1 / t \ ,„ r ^, k = l,...,z + l . (44) 

Cl,l (c(0)ct(0)e-' 3H /) o v ; 

Now from the Pauli principle we have the algebraic relation: 

c<{i)n(i) = (45) 



8 



which leads to the property 

c f (0) e- f3Hl = c f (0) . (46) 

Then Eq. (44]) takes the form 

Cl _ k _ (c(0)ct(0)K(0)]^-^ n 

Ci,i _ (c(0)ct(0)> ' I- 1 • V 



Now let us observe that Hq describes a system where the original lattice has been reduced to the central site (0) 
and to z sublattices, all disconnected among them and topologically equivalent to the starting one. Therefore, in the 
Ho-representation the correlation functions which connect sites belonging to disconnected graphs can be decoupled: 

</ {n (0)} g {n (p)})„ = (/ {n (0)}) (g {n (p)}) ( 
(g {n (p)} h {n (q)}) = (g {n (p)}) (ft {n ( 9 )}) ' ^ 

Here / {n (0)}, g {n (p)} and h {n (q)}, with p and q belonging to different sublattices, are any functions of the particle 
density. By means of such a property, Eq. (|47[) can be cast in the following form 

gi = < /[n"(0)] fe - 1 ) o , fc=l,...,*+l. (49) 

In Appendix C we show that 

([n a (0)] k } =F^[X] k=l,...,z (50) 
where F {z ^ [A] is a polynomial of order k in the variable X, defined as 

* = <« Q (0))o-^, (51) 



whose explicit expression is 



F (,,fc) [x] = E ^ =i ^ = ^ (52) 

here 6^ are some numerical coefficients defined in Appendix C. The previous analysis shows that all the properties 
of the system can be expressed in terms of only one parameter, A, defined by Eq. (|5ip . In order to determine this 
parameter we use the self-consistent equation 

C m = 1-k (1) , (53) 

where we required the translational invariant condition (n° (i)) = (n(i)). In order to exploit this equation we note 
that from (|3"9")l . by using the definition P7|) . we get: 

z+l 

lab = Ca.b + ^ ' ^am [fimp] 6 ^ Em Cp t0 . (54) 
m,p— 1 

Now, writing such equation for In and l\i and recalling that In = 1, lyi — K^- 1 ', we obtain 

z+l 

1 = Ci,i+ ^mPmp^e-^Cp,! (55) 

m,p— 1 
z + l 

= Ci )2 + ^ fiim[n m p] _1 e- fiEm C Pt2 . (56) 

m,p— 1 

Putting such expressions in the right hand side of the self-consistent equation ([53")) we get 

z+l 

Ci tl = Ci,i - Ci, 2 + ^ H lm [n mp ]- x e"^ (C p ,i - C p , 3 ) , (57) 

m,p— 1 
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which, by using Eqs. (|49)l and (|5"T)) . can be rewritten as 



z+l 



1 = (1 - X) + e- 0Em W {z ' m \ (58) 

m=l 

where 

z+l 

w ( Z ,m) = j2 n lm [n^r 1 (< [n a (O)]^ 1 > - < [n a (0)} p >„) . (59) 

The result (1ST))) allows us to express the function W^ z,m ^ in terms of the parameter X; indeed it can be shown that 

W^ = ^-_\)x m -Hl-Xr +1 - m m = W (60) 

yy(z,z+l) _ q 

Then, it is possible to write Eq. (|58]) as follows: 

' z- 1 



1 = (1-X) 



i+E e_w " u-i p m ~ 1(1 - x) 



(61) 



Recalling that E m = —\x — 4 (m — 1) J, /i = 2 (/i — z J) and making some algebraic manipulations such equation takes 
finally the form: 

X = (l-X)e^[l+ (e 4/3J - 1) X] . (62) 

Eq. (|62p is the main result of this Section; it allows us to determine the parameter X in terms of the external 
parameters /i, T, J (or h, T, J). 

Let us notice that, if we define a parameter x as 



e 2/3J 



l + ( e W_i)x' 

then Eq. (|62|) can be cast in the form 



(63) 



( 2/3,7 _ . 

f flPh = x *-i \e x) (g4) 



(xe 2 / 3 ' 7 - 1) 

which coincides with the one given by Baxter (see Ref. [25j, p. 53). In our case the quantity X has a definite physical 
meaning, it is the particle density of the first shell in the i?o- r ep re sentation. 

V. LOCAL CORRELATION FUNCTIONS AND RELATED PHYSICAL QUANTITIES 

The aim of this Section is to compute all the local correlation functions by expressing them in terms of the parameter 
X introduced in Eq. (|51[) . The calculation of the relevant physical quantities, that is particle density, magnetization, 
susceptibility, internal energy, specific heat and entropy per site, then easily follows. 

Let us start by recalling the results (|4"9"|) , ([50)) and (|5li)) which allow us to write the correlation functions in terms 
of the parameter X as follows 

Ci,k = C lt iYZ=i<>v' k ~ 1)xp > k = 2,...,z + l. (65) 
where C\ t \, due to Eq. (|46)) , can be expressed as: 

_ <c(0)ct(0)) o 

Cl4 " (e-^) ■ (66) 
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In order to compute C\ t \ let us observe that in the _ffo-representation c(0) satishes the equation of motion 



Then it is immediate to see that: 



*JU(0) = - M c(0). 



(67) 



(68) 



In order to evaluate the quantity ^ Hl / let us observe that, by means of the algebraic property [n(i)] m = n(i), 
we can write 



e -0Hj = e ^Jn(0) Zl =1 »(P) = f[ [1 + An (0) n (p)] 



(69) 



P =i 



where A = e 4/3J — 1. By using the property (|4"B")) and by recalling that in the Bethe lattice all sites are equivalent and 
the parameter X satisfies Eq. (|6"2")l , straightforward calculations show that: 



1 i a x 2 

( e ^ Hl )o = 1 + 1(1 + W 1] (n (0)) = {1 _x HePll + 1 y 



Putting dSHl and (J!| into (J66J) we finally get 



Ci. 



1 - X 



1 1 + AX 2 ' 



Now we are ready to calculate the particle density 



n=(n(0)) =1-C M = 



X (1 + AX) 

1 + ylA 2 ' 



the magnetization 



and all the correlation functions 



= <5(0)) = 2(n(0))-l = 



X (2 + AX) - 1 
1 + AX 2 



Ci,fe 



Y,Li a P ' XP > k = 2,...,z + l 



l+AX 2 ^p= 

Let us now switch to the calculation of the correlation functions 



K (k) _ /[„«» ( Q )]fe\ 

A( fc ) = (n(0) K (0)]* 



j.\ , A: — 1, z. 



According to the scheme given in Appendix D we have 

~ l + AX 2 ^ a P X 



P =i 



(l-X) + 



X{l+A) p 
(1 + AA) P_1 



1 y, fefc) XP +1 (l + Af 

1 + AA 2 ^ V (l + AXf- 1 ' 



The susceptibility per site can be calculated by means of Eq. ([73]) and has the expression: 

_dm _ !3{l~m 2 ) 



dh 



where we introduced the parameter p, defined as 



i-0 - i)p 



AX(l-X) 
1 + AX ' 



(70) 



(71) 



(72) 



(73) 



(74) 



(75) 



(76) 
(77) 

(78) 
(79) 
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This expression coincides with the one given in Refs. [26[ with p playing the role of the ratio 7^ of the eigenvalues 
of the second order transfer matrix V. 

Recalling the Hamiltonian ([1]), we obtain for the internal energy per site 

E(T) = ±(H) = -J [m 2 (l-p)+p]-hm (80) 
where we used the fact that the total number of points in the graph is [25[ 

N = 1 + z {z - l)^ 1 = l + z i - Z ~ z ^ ~ 1 . (81) 

q=l Z 

Once E (T) is known, we can directly calculate the specific heat, the free energy and the entropy (per site) by means 
of the formulas: 



f T E[T) — E (T*) 

F(T) = E(T*)-T V ' „ dT, (83) 

Jt* T 2 

where the value of E is given by l|80p and the limit T* — + is understood. 

VI. NON LOCAL CORRELATION FUNCTIONS AND RELATED PHYSICAL QUANTITIES 

In this Section we will calculate the relevant non local correlation functions; then we focus on the spin-spin one and 
on the related correlation length which we compare with the results existing in the literature [26j. We will show how 
our procedure allows us to evaluate also higher order non local functions with respect to the one given in Refs. [26| . 
Further technical details are presented in Appendix E. 

A. Two-point correlation functions 

Let us start by defining the correlation functions 

K^{j) = (K(0)] fe n(j)) (85) 
A«(j) = (n(0)[n a (0)]*n(i)) (86) 

where j is a site at a distance of j steps from the central site. Let us make for simplicity the choice that j belongs to 
the z-th subtree (but any subtree can be chosen) and let us focus first on the two functions: 

K^(J) = (n a (0)n(j)) (87) 
A(°)(j) = (n(O)n(j)). (88) 

We see that A^ ) (j) is a two-point correlation function which connects two sites which are j steps apart. Observing 
that (n (2) n (j)) and (n (p) n (j)) , z connect two sites which are j — 1 and j + 1 steps apart, respectively, it is 

immediate to see that the two correlation functions K^> (j) and (j) are related through the following relation 

K (D (j) = I A (0) {] _ 1} + i^l A (o) {j + 1} {89) 

z z 
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Let us now study the function A^ ) (J); it is immediate to see that 



A(°) (0) = (n(0)) =n 
A(°) (1) = (n (0) n (1)) = At 1 ) = n 2 + n (1 - n)p (90) 
A(°) (2) = (n (0) n (2)) = (zk< 2 ) - n) = n 2 + n (1 - n) p 2 



where we used the results of Appendices A and D [cfr. Eqs. (|D10[) - (|D12[) ] and we noticed that by means of (|D10[) the 



parameter p, defined by Eq. (|79p . can be expressed as p — X n(i-n) . On the other hand, in Appendix E we prove the 



following recursion relation 



G (j + l)-pG (j) = 1 [G (j) - P G (j - 1)] (91) 
p{z-l) 



where we defined 

G(i)=A(°)(j)-n 2 . (92) 
Then, the two-point density correlation functions A^ ) (j) for any j take the expression 

A(o) (a\ _ „2 

A(°)(j) = n 2 +n(l-n)p> ^ U ' ~ = p> . (93) 

n (1 — n) 

We are now in the position to calculate higher order correlation functions. By putting (|9"3"|) into (fBU)) we get 

K« (j) = n 2 + -n (1 - n) [p^ 1 + (z - 1) p> +1 ] , (94) 

z 

while, by putting the result <[93|) into Eq. (1E13|) of Appendix E we obtain: 

A (1) (j) = nA (1) + (1 ~ - {np^ 1 + [1 + n (z - 2)] p J + (1 - n) (z - 1) } . (95) 
By noting that the parameter p can be written as 

« * -i ("° (0))o , Qfi x 
n (n (0)) 

we see that it is always p < 1. Then the correlation functions A^ ) (j), A^ 1 ' (j), A^ (j) satisfy the ergodic theorem: 

lim^ooAC) (i) = (n(0)) (n(j))=n 2 

lim^oo A-W (j) = (n Q (0)> (n (j)) = n 2 . (97) 
lim^oo A« (j) = (n (0) (0)) (n (j)) = rcA« 

Recalling ([4]), we can evaluate from (|93p the spin-spin correlation function 

(S (0) S (j)) = m 2 + (1 - m 2 ) p> . (98) 

This expression coincides with the result of Refs. [26[. Now, by defining the correlation function: 

G s (j) = {S(0)S(j))-{S(0))(S(j)) (99) 

we obtain from (|98|) 

G s (j) = (1 - m 2 ) e-t (100) 

where the correlation length is defined as 



In 



(101) 



13 



B. Three-point correlation functions 

Following the same line of reasoning which led us to the two-point correlation functions, let us now calculate 
three-point correlation functions. Let us define the general three-point correlator as: 

T( fe ) (j,w) = (0) ct (0) [n a (0)f _1 n (j) n (w)\ = M^- 1 ) (j» - N^ 1 ^ {j,w) (k > 1) (102) 

where we introduce the new correlation functions 

M<*>0» - (K(0)] fc n(i)nH) (103) 

JV<*> (j,w) = (n(0)[n a (0)} k n{j)n(w)y (104) 

By j and w we denote two sites at a distance of j and w steps, respectively, with respect to the central site (0). Let 
us distinguish the two following cases: 1) j and w belong to the same subtree; 2) j and w belong to different subtrees. 

Case 1 

j and w belong to the same subtree, which we take as the ^-subtree, but any subtree can be chosen. By means of 
(jl2|, (J45J) and ((48)) and by noting that: 

( e -l3H l)o =Ci,i = l-n, (105) 

we can express T^ k > (j,w) as: 

T (fc) (j,w) = M^- 1 ' (j» - A^" 1 ) (j» = (1 - n) ^[n Q (O)f- 1 n (j) n (uj)^ . (106) 

Let us concentrate the attention on the two functions (j, id) = (n (0) n (j) n (w)) and (j, id) = 

(n a (0) n(j)n {w)). At first, we note that these functions are related through the following relation 

M« (j, w) = -NW (j - l, w - l) + Zzltf<P) u + i, w + i) . (107) 

z z 

Next, let us study the function (j, w); by recalling the definitions (|75|) it is immediate to see that 

(0, 0) = (n (0) n (0) n (0)) = (n (0)) = n 
N<® (0, 1) = AT(°) (1, 0) = iV(°) (1, 1) = (n (0) n (1)> = A (1) 
AT(°) (0, 2) = JV<°) (2, 0) = iV(°) (2, 2) = (n (0) n (2)) = -X^ (zk' 2 ' - n) ' (108 ^ 
JVC°) (1, 2) = AT(°) (2, 1) = (n (0) n (1) n (2)) = jlj (zA( 2 ) - A (1) ) 

Recalling the relations (|D12[) in Appendix D and the expression p = ^^Z") i the correlation function A^ -* (j, w) can 
be written in the closed form 

N ( °) (j :W ) = n 3 + n 2 (1 - n) (p> + p w ~ j ) + n (1 - n) 2 p w , j,w = 0,1,2. (109) 

In order to evaluate (j,w) for all values of j and w we need a recursion formula as the one in Eq. (I91|) . which we 
now derive following the same steps outlined in Appendix E for the function A^ -* (j). 
In the i?o- re P resen tation (j,w) can be written as 

iVW (j, w) = („ (0) n (j) n («;)) = V ^. (110) 

Recalling that (cfr. ([Ml)) e~^ H ' = J]p = i I 1 + An (0) n (p)], by making use of Eqs. flgSJ) and ((701), as well as of the 
equation (|62p for the parameter X , we have 

^ (0) U, w) = 1+ X AX2 (n (j) n H (1 + An (*))>„ . (Ill) 
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Let us now calculate the function (J, w) = (n (0) n a (0) n (j) n (w)). By following the same procedure, we obtain 
for j, w > 1 : 

iV«(j» = {Z = 1} j, ; ^± ^ ^ (n (j) n M (1 + An (z))) 



+ - 



z (1 + AX) (1 + AX 2 ) 
11(1 + 4) 



(n (z)n(j)n(w)) . 



(112) 



z 1 + AA 2 

For simplicity, let us restrict the analysis to the case where w follows j (i. e. we can write w = j + r, where r is the 
number of steps necessary to go from w to j). In this case we have Af(°) (j, to) = (r) and Eq. (flQgj) for fc = 1 
gives 

AC) (r) - AT(°) (j,j + r ) = (1 - n) (n (j) n (j + r)) . (113) 

By combining ()113|1 and (jllll) we can express the two quantities (n (j) n (j + r)) and (n (j) n (j + r) n (z)) in terms 
of JV(°) (j,j + r) and A^ ) (r) as follows: 

1 



(n{j)n(j + r)) Q = 
(n(j)n(j + r)n(z)) = 



(1-n) 
l + AA 2 



A (°) (r)-J\T<°> (jJ + r) 

1 1 



AA (1 - X) 

By using the relation p06[) for fc = 2 and by observing that 



A (0) (j,i+r)- 



A(l-n) 



A<°> ( r ) . 



1 z — 1 
(n a (0) n (j) n (j + r)) = - (n (j) n (j + r) n (z)) H X (n (j) n (j + r)) . 



z z 
we can express the function M^- 1 ' + r) in terms of A^ ) (jj -f r ) and A^°) (r) as 



(114) 



(115) 



M (1 » (jj+r) 



(*-l)p , 1 
z zp 



N {0) (j,J+r) + 



(z-l)p 1 

z zp 



,(l-p)A(°) ( r ). 



Recalling that [see (p?3"l) ] A^°) (r) = ri 2 + n (1 — n)p r we obtain from (|107[) the relevant recursion rule 

1 



Ar ( °) (j + i, w + i) _ p AT ( °) (j, to) = 



AT ( °) (j, to) - pN {0) {j - 1, to - 1) 



(116) 

(117) 
(118) 

(119) 
(120) 

(121) 

this expression agrees with the one given in Ref. [33j, where it was calculated for the one dimensional case (i. e. 
z = 2). It is interesting to notice that the expressions of the spin correlation functions [cfr. (f9"8"| and (|121[) ] depend 
on the coordination number z only through the parameters p and m. 

We are now in position to calculate the correlation functions AfW (j, w) and A^ 1 ) (j, to). Straightforward calculations 
give 

M« (j, to) = n 3 + -n (1 - ri) {n (p>~ 1 + p 10 ^') + (1 - ri) p^ 1 + («-!) [n (p 3+1 + p w - j ) + (1 - n)p w+1 ] } , (122) 



p(z-l) 

where we defined 

Ar(0) (-,; w ) = jvCD) (j, „,) - n 3 _ n 2 ^ _ n ^ pW -j 

By recalling the result (|109|) , it is easy to see that the following expressions hold for any j and to 

AT(°) (j, to) 
n (1 — ri) 
A ( °) (j, w ) = n" 

By using the transformation (j4|) we obtain from (|120j) the expression of the three-spin correlation function: 

(S (0) S (j) S (to)) = m 3 + to (1 - to 2 ) (y - p w + p 1 "^') ; 



= np> + (1 - n)p™ 

^ 3 + n 2 (1 - n) {p> + p<"- J ) + n (1 - n) 2 p u 



(j, w) = n 2 A (1) + n 2 (1 - n) [np 10 ^' + (1 - n) + 



nil — ri) 



i 2 p j ~ 1 + (1 - rt) 2 ^ + n (1 - n) {p w ^ + p 7 ') 



(123) 



n (1 — ri) (z — 1) 



n 2 p> + (1 - n) 2 p"' +1 +n(l - n) + p w ) 



Let us notice that all the correlators A^ ' (j, to), AfW (j, to) and A^^ 1 ^ (j, to) satisfy the ergodic theorem. 
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Case 2 

Let us now switch to the second case: j and w belong to different subtrees, which we take as the z- and (z — 1)- 
subtree, respectively. By performing the same steps which led to Eq. (|10€>|) and by noting that in such a case 
M(°) (j,w) = A(°) (j + w) we obtain the relation 

N<® (j,w) = A<°> (j +w)-(l-n) (n (j)> (n H) . (124) 
Recalling the expressions (fETTj) for (n (i)) and for A< 0) (j + to), we obtain for any j and w 

N {0) (j, w)=n 3 + n 2 (1 - n) {p> + p w ) + n (1 - n)V" Hu - (125) 
To calculate higher order correlation functions we observe that 

A/« (j, w) = -N® (j - 1, w + 1) + i^ (0) (j + 1, u; - 1) + ^^AT(°) (j + l, w + l) . (126) 
Putting ([125]) into ([126]) we have 

M (1) (j, u>) = n 3 + -n (1 - n) [n (p 3 '" 1 + p w ^) + 2 (1 - n) p ]+w + 

(z-l)n (p> +1 + p w+1 ) + (z - 2) (1 - n) p 3+tu+2 ] . (127) 

In order to calculate JVW (j, w) let us observe that 

7V« (j, w) = M « (j, w) - T< 2 > (j,™) = M« (j, w) - (1 - n) (n a (0) n (j) n («;))„ . (128) 
which, by noting that 

(n a (0) n (j) n («,))„ = ± (n (z) n (i)) (n (w)) + ± <„ (* - 1) „ (w)) (n (j)) Q + Z -^X (n (j)) (n (w)) , (129) 
and recalling (jElip and (jl2T[) . becomes 

N<V (j» = n 2 A« + [(z ~ 2)n+1] n 2 (1 - n) (p* + + ^^n 2 (1 - n) 2 (p> +1 + p w+1 ) 

z z ' 

+ n3(1 ~ n) (p*" 1 : ) + -n 2 (1 - n) 2 ^'+— 1 + { ^—^n 3 (1 - n) p> 

z v ' 2 z 

t [2(l-2n)-(* -4)n 2 ] (z-2) , 2j+M+2 

-7i (1 — n)p J H u(l — n) p J (1«j0) 



z 

for any j and it). Also in this case the ergodic theorem is satisfied. 

VII. RESULTS 

In the previous Sections we have shown that all the properties of the system are expressed in terms of the correlator 
X = (n a (0)) Q . This quantity is determined in terms of the external parameters J, T, h by solving the equation {62}, 
which is a polynomial of order z in the variable X . In this Section we discuss the solutions of Eq. (|62[) and present 
the results obtained for various properties: the magnetization, the susceptibility, the specific heat, the free energy 
and the entropy. We shall discuss separately the cases of zero and finite magnetic field, by focusing the analysis to a 
ferromagnetic coupling (i. e. J > 0). 

A. Zero magnetic field 

In the case of zero magnetic field we have /j, = —2zJ; then, it is useful to define K = e 2 ° J so that the equation (|6"2")) 
takes the form 

XK Z = (1 - X) [l + (K 2 - 1) X] 2-1 . (131) 
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It is easy to see that 



x TTk l^Tl (132) 



is always a solution of the equation (|131|) for any value of the coordination number z. By putting (|132p into l|72p and 
(173")) we have 

n = (n (0)) = |, m = (S (0)) = . (133) 

The particle density and the magnetization do not depend on the temperature and on the coordination number z. 
This is a manifestation of the particle-hole symmetry, when we recall the scaling law (|10[) for the chemical potential. 
But (|131[) may admit other solutions which break the particle-hole symmetry. In particular, let us study if there is a 
critical temperature T c such that the magnetization is different from zero for T < T c . In order to determine T c let us 
expand (|131|) in power series of X around the point given by (|132|) . At first order we obtain 

X- Y ^\[z(K-l)-2K]=Q. (134) 

Therefore, besides the solution (| 1 32[1 there may be other solutions when K = . Such an equation shows that there 
is a critical temperature T c , given by 



2J = k B T c \n\-- -j , (135) 

such that for T < T c we may have solutions which spontaneously break the particle-hole symmetry and exhibit a 
magnetization different from zero. Let us notice that the case z — 2 (i.e. the one-dimensional chain) gives T c = 0. Let 
us also point out that Eq. (|135p admits a solution only when J > 0. For negative J (i.e. antiferromagnetic coupling) 
there is no solution. If z is even, the equation (] 1 3 1 [) admits another solution 

which also gives the results quoted in (|133|) . However, such a solution describes an unstable system: the energy is a 
decreasing function of temperature and the parameter p is larger than one. This solution will be disregarded in the 
following. 

Generally, for z > 2 we have the following situation: 



• z even 



2 solutions corresponding to n = i and m = 
T < T c ^ 2 solutions corresponding to ± m ^ 

the remaining roots are complex 
T > T \ 2 solutions corresponding to n — g and m = 
1 the remaining roots are complex 



• z odd 



1 solution corresponding to n — \ and m — 
T < T c ^ 2 solutions corresponding to ± m ^ 

the remaining roots are complex 
j, j, . 1 solution corresponding to n = 5 and m = 
1 the remaining roots are complex 

By considering the following items: (i) the broken symmetry solution (i. e., m 7^ 0) has a free energy lower than 
the one corresponding to the symmetric solution, (ii) the solution (| 136[) is disregarded because not physical, (iii) the 
two solutions corresponding to ±m are physically equivalent, (iv) all the complex solutions are disregarded, we can 
assert that the equation (|13ip admits only one solution of physical interest. 

For T > T c we have the following results: 

X = i [1 - tanh (j3J)} n = | = | [1 + tanh (/3 J)] p = tanh {(3 J) 

m = £ = -Jtanh(/3J) C = k B [/?Jsech (/3J)] 2 X = ^^tl^j) ' { ' 
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For T < T c the breaking symmetry solution depends on z. We shall present results for z = 3 and z — 4. For the 

' sTa _ _2_ 
J In 3 



case of z = 3 the critical temperature is given by kB J a — r^r ~ 1-82048. The solution of (|131[) is 



_ (A-+l)(A'-2)+AV(K+l)(K-3) 



(A-+l)(g-2)+A-y'(A+l)(A'-3) _ if ^/(jf +l)(/f-3) M38) 

~ " ■ l){K-2) m - " 

1 , _ 40 K 



2(K+l)(K-2) 771 ~ (K+l)(K-2) 

y K-l A (if-3)(K-2) 2 (K+l) 



For the case of z = 4 the critical temperature is given by kB J c = ~ 2.88539. The solution of (| 1 3 1 [) is: 

Y _ K 2 -2+KVK 2 -4 

A - 2(g*- l) 

II — 2(K 2 ~2) " L ~ K 2 -2 ■ 

„ _ 1 v _ 

^ i?2=T A (K2- 4 )(A' 2 -2) ;! 

For all values of z the parameter A'- 1 -' and the internal energy E can be calculated by means of the expressions 

AW = n\p + n(l-p)], E = AJn (1 - n) (1 - p) - 2J . (140) 

In Fig. 1 we plot the magnetization per site m = (S (0)) as a function of the temperature, expressed in units of J, 
for the values of the coordination number z = 3 and z — 4. As expected, the magnetization decreases by increasing 
T and vanishes at the critical temperature T c , determined by (|135p . By expanding the parameter K = e 2/3J around 
the critical temperature T c : 

K = K c [l + bt] + (t 2 ) (141) 

where 

t = ^, if c = ^, & = ln(-0 , (142) 



T c ' " c z-2 > " \ 2-2 

we can easily show that close to T c the magnetization behaves as 



^+0 if z = 3 
m = { %- > A (143) 



with critical exponents /? = i, in agreement with Refs. [25|, [26[. The behaviour of the parameter p = 1 — ^ ■ 



2> 111 o-Biccmciiu vviuu hud. j^iijjLsijj. me uraianuiu ui unc paiamcuu ^ — a. (n(0)) 

as a function of % is shown in Fig. 2. By increasing the temperature, p first increases up to the maximum value 
{p)t=t = ~ ^T' ^ nen decreases. It is always p < 1: this condition implies the ergodic behaviour of the spin correlation 
functions, when we recall the results of Section 6. We notice that for T > T c the value of p is the same for all values 
of z. 

In Fig. 3 we report the temperature dependence of the spin susceptibility per site x — (^)h— n' This quantity 
diverges at T = T c with critical exponents 7 = 7' = — 1: 

{ (T e — T) 1 rri rp 
b(z-2) 1 > ±c 

To the contrary of p, above T c the susceptibility changes with z, as can be seen by Eq. (|75|) . The specific heat per 
site is reported in Fig. 4 as a function of the temperature. We observe a jump in correspondence of T c , with critical 
exponents a — 0, as expected for a second order phase transition. It can be shown that the jump AC at T = T c is 
given by: 

4J 2 K C 




(1+K C ) 2 T 2 ~ 8 [ ln3 l Z ~ 3 

4.7 2 jf e 
{l+K a ) 2 T 2 



AC -< *n.J 2 u 2 K, _ 2[ln2] 2 z = 4 • (145) 



We notice that the jump decreases with z, and that above T c the behaviour of the specific heat does not depend on 
the value of z. 
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The temperature dependence of the internal energy E, of the free energy F and of the entropy S is shown in Figs. 
5, 6 and 7, respectively. We observe the different behaviour at T c : F is a smooth function, while E and S exhibit a 
drastic change. This behaviour shows that at T c we have a second-order phase transition. Also, we note that above 
T c the internal energy does not depend on z, while the free energy and the entropy depend on z. 

In Fig. 8 we plot the spin correlation function A^ ^ (j) — (n (0) n (j)) versus the distance j for z = 3 and several values 
of the temperature, chosen below and above the critical temperature. We clearly see that a long-range ferromagnetic 
order is established below T c . 

B. Finite magnetic field 

For finite magnetic field the equation (|62p for the parameter X can be written as 

(1- X)H[l+ (K 2 -1) X]^ 1 - XK Z = (146) 
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FIG. 4: The specific heat C is plotted against T/J for z = 3 and z = 4 and zero magnetic field. 



where we put K — e 2/3J and H — e 2>3h . For H ^ 1 this equation does not admit a general solution for any value of 
the coordination number, and we must discuss case by case. For z — 2 the solution is 



X 



K 2 (H - 1) - 2H + K^jK 2 + H 2 K 2 - 2H (K 2 - 2) 
2H (K 2 -l) 



(147) 



which describes the well known solution of the one-dimensional spin-^ Ising model. The other root of Eq. (|146p 
corresponds to a physically unstable system and is disregarded. 

For z — 3 it is possible to show that there is a critical temperature T c (h), depending on the magnetic field, such 
that for T < T c (h) there are three real and unequal roots, while for T > T c (h) there is one real root and two conjugate 
imaginary roots. T c (h) is determined by the following equation: 



8K 3 H = K A 



18K 2 - 27 - (K 2 - 9) 3/2 VK 2 - 1 



(148) 



For z = 4 there are four real and unequal roots for T < T c (h), while for T > T c (h) there are two real unequal roots 
and two conjugate imaginary roots. 

In Fig. 9 we report T c (h) as a function of the magnetic field for z = 3,4. By increasing \h\, T c (h) decreases from 



the value fcsT c = 



2.7 



Ma) 



[cfr. Eq. (|135[) ] and vanishes at \h\ = J (z — 2) 
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FIG. 9: The temperature T c (h) is plotted against the magnetic field h for z = 3,4. 



For z = 3, below T c (h) the three real solutions have the following behavior: one solution satisfies the particle-hole 
symmetry (i. e. (n) = i at )U = —2zJ) but corresponds to an unstable system (for example the compressibility is 
negative). The other two solutions violate the particle-hole symmetry and describe a finite magnetization for any 
value of the magnetic field, in a direction parallel and anti-parallel to h, respectively. By disregarding the unstable 
solution and by picking up the one which describes a magnetization in the direction of the magnetic field, we can 
assert that in the entire plane (h, T) Eq. (|146[) admits only one physical solution given by 

X = cos (0) - § T <T C (h) 



X=^/f = ^+^/f = ^-l T>T c (h) (149) 

where 

_ K 3 (KH-3) , = K 3 (2HK 3 -9K 2 +27) 

a 3 ff(tf'-i) a 27H ^ K )-\l N . (150) 

r - 3-K 2 a _ 1 „ os -l ( t!^l \ K ' 

C -TCZ—1 t7 - 3 COS ^ 2a 3/2^ 

Similar situation holds for z = 4. Below T c (ft,), two solutions satisfy the particle- hole symmetry but correspond to 
an unstable system; the other two solutions violate the particle-hole symmetry and describe a finite magnetization, 
parallel and anti-parallel to h, respectively. Above T c (h), among the two solutions, only one has a physical meaning. 
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Once X is known, we can calculate the various properties by using the formulas given in Sections 5 and 6. The 
behaviour of the magnetization, the susceptibility and the specific heat as functions of the temperature is reported 
in Figs. 10, 12 and 13, respectively, for several values of the magnetic field. At low temperatures the system is fully 
polarized by any finite magnetic field. By increasing temperature, the magnetization decreases and tends to zero in 
the limit T — > oo. The different behavior of the magnetization, below and above the critical temperature T c , is shown 
in Fig. 11, where m is plotted as a function of the magnetic field for several values of T. For finite h the susceptibility 
and the specific heat do not exhibit a discontinuity, there is a peak at a certain temperature T* which increases with 
h. 

1.2 
1.0 
0.8 
0.6 
0.4 
0.2 



m 




h=0 

h=0.1 

h=0.5 

h=l .0 



0.0 \- z=3 
-0.2 



J I I I I I I I I I I I I I I I I I I I I I I L 



1 2 3 4 5 

T/J 

FIG. 10: The magnetization m is plotted against T/J for z = 3 and several values of the magnetic field. 




0.5 1 1.5 2 

h/J 

FIG. 11: The magnetization m is plotted against h/J for z = 3 and several values of the temperature. 



VIII. CONCLUSIONS 



In this paper we have studied the ferromagnetic spin-i Ising model in the presence of an external magnetic field 
on the Bethe lattice by means of the equations of motion method and of the Green's function formalism. First, 
it has been shown that there exists a mapping between such a model and a fermionic one built of one species of 
particles localized on the sites of a Bethe lattice and interacting via an intersite Coulomb interaction. Then, an exact 
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FIG. 12: The susceptibility \ is plotted against T/J for z = 3 and several values of the magnetic field. 
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FIG. 13: The specific heat C is plotted against T/J for z = 3 and several values of the magnetic field. 

solution has been found for the Bethe lattice with any coordination number z in terms of a finite and complete set 
of eigenoperators and eigenvalues of the fcrmionic Hamiltonian. This solution allowed us to write exact expressions 
for the corresponding Green's function and correlation functions which depend on a finite set of parameters to be 
determined in a self-consistent way. Such parameters have been exactly fixed by means of algebra constraints. So 
local and non local correlation functions have been calculated up to four point ones, together with the corresponding 
physical quantities, i.e. the particle density, the magnetization per site, the susceptibility, the correlation length, the 
internal energy per site, the specific heat and the entropy. All the results are in perfect agreement with the ones 
existing in literature. The results are discussed in great detail for the first two values z = 3,4 of the coordination 
number of the Bethe lattice with and without external magnetic field. Our procedure allows also to generalize the 
known results, as it has been explicitly shown in the case of non local correlation functions, where new expressions, 
no previously reported, for the three-point spin-spin correlation function have been obtained, together with a general 
calculation scheme which can give rise to higher order correlators. 
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APPENDIX A: CALCULATION OF [n a ] k 
Let us recall the following definition 

n a 00 = - (m +n 2 + ■■■ +n z ) (Al) 

where n p (p = 1, z) are the first nearest neighbors of the site %. Then we have to calculate the power 

[n a (i)] k = ^(n 1+ n 2 + ... + n z ) k . (A2) 

By considering the algebraic property [n p ] m = n p (to = 1,2,...), after some easy but lengthy calculations, it can be 
shown that 



[» a w] fc = 4 EM*) (A3) 
z p=i 

(z) 

where Np ' are the operators 

z 

N { p z) = E n h n h ...n lp (A4) 
Zi<i 2 <.--<i P =l 

and are some positive integer numbers, defined as 



b[ k) = 1 



(A5) 



(A6) 



,(fe) _ v^/l-l r l„-3-l r lp-!-l f k \ f ll \ f l2 \ ( l P- 3 \ ( l P- 2 \ <rx ^ 11 

o P - 2^ 1=P -i 2^ a=P -2 U 3 = P -3 - 2^ p _ 2 =2 2^ p _ 1= i ^ ii i 2 J V J'" V ^-2 / ^ J (P ' 

(k) 

It is important to notice that b p — for p > k. In particular for the first values of p we get: 

4 n) = 2 (2"- 1 - 1) 
4 n) = 3 (3™- 1 - 2™ + l) 

= 4 ( 4 n-l _ 3 n + 3 . 2 n-l _ l) 

6< n) = 5 (5™- 1 - 4™ + 2 • 3" - 2™ +1 + 1) 
The above results are valid for any lattice with coordination number z. 

APPENDIX B: CALCULATION OF THE COEFFICIENTS A$ 

Given the results in Appendix A, we can write 

[n t "(i)] fc =E^ ) [« a Wr (Bl) 

m— 1 

where the coefficients A^m are some rational numbers which must satisfy the relations 

V 2 A( k ) — 1 

4 «rr^-i ^ (B2) 

The first relation follows by putting n a (i) = 1 while the second can be derived by considering the case 1 < k < z. 

(k) 

Indeed we must calculate the coefficients Am only for k > z + 1 and m = 1, z. By noting that for k > z + 2 we 
can write 

E 4? K HT = E A ^ 1] K (0] m+1 - (B3) 

m — 1 m— 1 
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(B4) 



the following recursion rule can be established 

At ] = + Af- 1] At +1) (m=l, ...,*) At 1] =0- 

This rule implies that we must calculate only the z coefficients Am +1 \ (m = 1, z), which can be done by means of 
the relation (IBID evaluated for k = z + 1: 



K«r +1 = E^ +1) [» Q (or. 

m— 1 

By using the results given in Appendix A to rewrite (|B5|) . we obtain the equation 



fc=i 



fc=i 



which, by noting that the operators Nf, are linearly independent, takes the form 

E™ =fe At +1) z^ m b^ - = (k = 1, z) ■ 

Such equations give rise to the iterative solution 

,(2 + 1) _ b^_ +1) 



A 



(2+1) i_ 

*" 2 ~ Z 3 b (- 



(2+1) _ 1 



26?* 

,(2 + 1)" ,(2 + 1) AZ) 

°z-l ~ A * ZO z-l 



b^-A^zHl^-A^zb^ 



(B5) 



(B6) 



(B7) 



(B8) 



(2 + 1) _ 1 



I- A, 



(2+1) ,2-1 _ A (z+\) vZ -2 



/((2+1) 3 4(2 + 1)^2 4(2 + 1). 

/1 2-2 Z _ ^2-1 Z _ ^2 



As an example, we give the values of the coefficients A m for the first values of z 



z 


= 2 : 


4 3 > = 


1 
2 

2 
9 


4(3) _ 3 
^2^2 














z 


= 3 : 




,(4) _ 11 

vr 9 




2 


Af- 








z 


= 4 : 




3 

32 


4(5) 25 
A 2 — 32 


4 5) = - 


35 

" 16 


5 
2 


4(6) 




z 


= 5 : 




24 
625 


4(6) _ 274 
^2 — 625 


4 6) = 


9 

5 


4 6) = 


17 

5 


= 3 


z 


= 6 : 


4 7 > = 


5 
324 


4(7) _ 49 
^2 — 216 


4 7) = - 


203 
162 


4 7) = 


245 
72 


4 7) = 


175 
36 



<(7) _ 7 
4 6 —2 



(B9) 



APPENDIX C: CALCULATION OF ([n a (0)] k ) Q 
By means of the results given in Appendix A we can write that: 

(i» a ®r)r?£ b P( N '% (C1) 
p=i 

The particular topology of the Bethe lattice allows us to decouple the correlation functions in the ffo- re P r esentation 
(see Eq. (|4"g)) ) and to obtain 

(K Z) ) = E (^n h ...n lp ) o = [(n a (0)) o } p E = (pW(°)>oF. (C2) 

ll<l 2 <---<l p = l 1 1 <1 2 <---<1 P = 1 ^ ' 

It follows that 

(K(o)] fe > --iE(;)^ p ( c3 ) 

p— 1 ^ ^ 

where X = (n a (0)) (see Eq. ([5T])). 
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APPENDIX D: CALCULATION OF A (fc) AND k w 

In this Appendix we will calculate the local correlators X^ and K^ k ' in terms of the parameter X defined in Eq. 
(ISTj) . We recall the definitions given in Section 5: 



A( fe > = (n(0) [n a (0)]' 



Let us start with X^ which, according to f4"2"| . can be written as 



'n(0){n«(0)] k e^ H ') 



By recalling Eq. (|69|) and the recursion rule (|A3|) . we have 



n(0) [n a (0)] e 



o z 



i=l 



Furthermore, by recalling Eq. (|A4[) it is immediate to show that 



(1 + Af X p (1 + AX) 



z-p 



(N^Y[[1 + A ni }) = 

\ i=l / 

By putting ([D3]) and (jD4f into (|D2jl and recalling the results ([68]) and (J70|) we obtain for k > 1 



A (fc) = _0__^_ e /3M 



1 + AX 2 z K ^ \ p 

P =i 



sZ-p 



Use of the equation (|62|) for the parameter A" allows us to rewrite (|D5j) under the form 



1 f- jf , fc ) ^ +1 (l + Af 



P =i 



where 



(*,fc) _ J_ ft (fc) f * \ 



In order to calculate , let us start from the equation: 

^=C Ii+ i+A (fc) . 

By putting together (I74| and (|D6|) . we obtain: 



l + AX 2 ^ p 



P =i 



(1-X) + 



X(l + Af 
(l + AXf 1 



In particular, from (|D6|) and (|D9[) we have 

n _ X 2 (l + A) 



A (1 



1 + AA 2 



z z (1 + AX 2 ) (1 + AX) 



(Dl) 



(D2) 



(D3) 



(D4) 



(D5) 



(D6) 



(D7) 



(D8) 



(D9) 



(D10) 
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and 



X (1 + AX) 



K 



(2) 



1 + AX 2 

-(i) 



z 



z 1 + AX 2 

In closing this Appendix, wc note the following useful relations 



(l-X) 



X(l + A) 2 
(1 + AY) 



where we used n = 



(2) 



1 Z-l 

—n H 

z z 



A' 2 



A (1 



(n 2 -A«)^ 
n (1 — n) 



lA , »-l (n-AW) s 
z (1 - n) 

l[A«] a 



(Dll) 



(D12) 



APPENDIX E: CALCULATION OF A (0) (j) AND (j) 

Let us recall the following definitions given in Section 6: 

K W(j) = ([n a (0)} k n(j) 



AW (i) = (n(0) [n Q m k n(j) 



(El) 



where j is a site at a distance of j (j > 2) steps from the central site. Let us make for simplicity the choice that j 
belongs to the z-th subtree (but any subtree can be chosen) and let us define the correlation function of the composite 
fields as 



D« (j) = (c (0) ct (0) K (O)f- 1 n (j)) (k > 1) 

By means of the commutation relations we note that 

£,(fc) (i) = K (k-i) {j) _ A (fc-D {j) t 

while, by using (|42]) and the algebraic relation (|46|) . we have 



£)(*) (7) = ^ L° 

By means of the properties of the correlation functions in the ifo-representation (see Eq. we get: 

, M (c(0)ct(0)) n / k , 

\ e /o 

= C M (K (O)]^ 1 [n (i)]) Q = (1 - n) (K (O)]^ 1 [n (j)]) Q > 
where we used Eq. (|53|) . By putting together (|E3|) and (|E5|) we obtain 

K (k-i) {j) _ A (k-i) (i) = ( i_ n) ( [n " (Q)]^ 1 [„ , 



which in particular for k — 1 reads 



(E2) 



(E3) 



(E4) 



(E5) 



(E6) 



(E7) 
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Let us notice that (n (j)) depends on the site j because the i?o-representation lacks of translational invariance. Let 
us now start to calculate the function A^ ) (j) — (n (0) n (j)), which in the iJo-representation can be written as 



A(°) (j) = (n (0) n (j)) = V 1 \ y> 



By recalling that (cfr. e~ 0Hl = n* = i [1 + An (0) rn] we have 



(n (0) n (j) e-^)o = ( n (°)>o C 1 + ^0 K» 0')>o + 4 0')> 



oJ • 



(E8) 



(E9) 



Putting (|E9|) into (|E8|) . recalling the results and ([TO)) , using the equation (jo^l) for the parameter A, we obtain 



A(°) (j) = 



A 



1 + AA 1 



[<»Cj)>o + ^<»« n O')>o] 



(E10) 



By combining (|E7|) and (|E10j) we can express the unknown correlation functions (n(j)) and (n z n{j))o in terms of 
the two-point correlation function A^ ^ (j) as follows 



(n(j)) 
( n zn(j)) 



1 



(!-») L- Aro «fl 
1 + AA 2 1 1 



AA 



A (1 - n) 



A(°) (i) 



1 n 
A (1 - n) ' 



(Ell) 



Let us now calculate the function Ad) (j) = (n(0)n™ (0)n(j)) by following the same procedure we adopted for the 
calculation of A^ ) (j). Recalling (JHHJ) an d ([70]) . and the basic equation ([52")) we get 



A« (j) 



(z-1) A 2 (l + A) 
z (1 + AA)(1 + AA 2 ; 
1 A (1 + A) 



[(« (i)) + A (^«(i))o] 



+ 



{n z n{j)) Q . 



z 1 + AX 2 

Now, by using the relations in (|E1 1|) it is possible to express the correlation function A^ 1 ) (j) in terms of A^ ) (j) 



(E12) 



A d) (A - ^(1 + A) 1 n (z-l)X(l + A) m 1 (1 + A) (0) 

Uj zl + AA 2 A(l-n) z (1 + AA) U) zA(l-X) [J) ' 

By using the relation (|E6|) for fe = 2 we are now in position to calculate the function K^> (j) as 

^ (1) (j)=A« (j) + (l-n) (n a (0)n(j))o. 

which, by observing that 

K (0) n (j)) = i (n z n (j)) + — A (n (j)> , 



and recalling (|E11[) and (|E13|) . becomes 



tf«(j) = - A(°)(j)-n 2 
zp 



(z-l)p 



A(°) (j) 



(E13) 



(E14) 



(E15) 



(E16) 



where we made use of the relations (fT2"j) and (1751) . 
Recalling now the recursion relation (cfr. (|89[)) 



(J) = -A(°) (i-l) + — AW y + 



and putting that together with (|E16|) . we finally obtain the relevant recurrence relation 

1 



G (j + l)- P G (j) = — — - [G (j) - P G (j - 1)] 

P{Z 1) 



(E17) 



(E18) 
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where we defined 

G (j) = A<°> (j) -n 2 = (n (0) n (j)) - n 2 . (E19) 
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